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ABSTRACT 

We present continued radio observations of the tidal disruption event Swift J164449. 3+573451 extending to 
(5f w 216 days after discovery. The data were obtained with the EVLA, AMI Large Array, CARMA, the SMA, 
and the VLBA+Effelsberg as part of a long-term program to monitor the expansion and energy scale of the 
relativistic outflow, and to trace the parsec-scale environment around a previously-dormant supermassive black 
hole (SMBH). The new observations reveal a significant change in the radio evolution starting at St I month, 
with a brightening at all frequencies that requires an increase in the energy by about an order of magnitude, 
and an overall density profile around the SMBH of p cx r~^^^ (0.1 - 1.2 pc) with a significant flattening at 
r « 0.4-0.6 pc. The increase in energy cannot be explained with continuous injection from an L cx f"^/^ tail, 
which is observed in the X-rays. Instead, we conclude that the relativistic jet was launched with a wide range of 
Lorentz factors, obeying E{> Tj) oc Ty^^. The similar ratio of duration to dynamical timescale for Sw 1644+57 
and GRBs suggests that this result may be applicable to GRB jets as well. The radial density profile may be 
indicative of Bondi accretion, with the inferred flattening at r 0.5 pc in good agreement with the Bondi radius 
for a ~ few x 10^ Mq black hole. The density at ^ 0.5 pc is about a factor of 30 times lower than inferred 
for the Milky Way galactic center, potentially due to a smaller number of mass-shedding massive stars. From 
our latest observations (St « 216 d) we find that the jet energy is £j iso w 5 x 10^-^ erg (Ej « 2.4 x 10^' erg for 
9j = 0.1), the radius is r « 1.2 pc, the Lorentz factor is Tj « 2.2, the ambient density is « « 0.2 cm"-', and the 
projected angular size is rproj ~ 25 /ias, below the resolution of the VLBA+Effelsberg. Assuming no future 
changes in the observed evolution and a final integrated total energy of Ej w 10^^ erg, we predict that the radio 
emission from Sw 1644+57 should be detectable with the EVLA for several decades, and will be resolvable 
with VLBI in a few years. 
Subject headings: 



1. INTRODUCTION 

The discovery of the unusual 7-ray/X-ray transient 
Swift J164449.3+573451 (hereafter, Sw 1644+57), which co- 
incided with the nucleus of an inactive galaxy at z = 0.354, has 
opened a new window into high-energy transient phenomena, 
with potential implications to our understanding of relativis- 
tic outflows in systems such as gamma-ray bursts (GRBs) and 
active galactic nuclei (AGN). The prevailing interpretation for 
this event is the tidal disruption of a star by a dormant su- 
permassi ve black hole (SMBH) with a mass of Mbh ^ 10^- 
. 10^ (iBloom et al."201 It lBurrows et al.ll20Tl [ Levan et ajj 
■ 1 201 U Z auderer et al. 20 0; but see iKrohk & P iran 2011]; 
' lOuved et al.ll201 It lOuataert & Kasen 201 1 for alternative ex- 
planations). The argument for a tidal disruption origin is 
based on: (i) a positional coincidence (< 0.2 kpc) with the 
host galaxy nucleus; (ii) rapid time variability in 7-rays and 
X-rays (< 10^ s), which requires a compact source of < 0.15 
AU, a few times the Schwarzschild radius of a ~' 10^ Mq 
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black hole; (iii) high 7-ray and X-ray luminosity of ^ 10 
erg s"', which exceeds the Eddington limit of a 10^ Mq 
black hole by 2-3 orders of magnitude; (iv) a lack of previous 
radio to 7-ray activity from this source to much deeper limits 
than the observed outburst, pointing to a rapid onset; and (v) 
long-term X-ray luminosity evolution following Lx oc f"^'''', as 
expected from the fallback of tidally disrupted material (e.g., 
lRees|[T988HStrubbe & Ouataertil2009i) . 

Equally important, Sw 1644+57 was accompanied by 
bright radio synchrotron emission, with an initial peak 
in the millimeter band (/v ~ 35 mJy) and a steep spec- 
tral slope at lower fre quencies indicative of self-absorption 
(IZauderer et al.l 1201 l l hereafter, ZBSll). The properties 
of the radio emission established the existence of a rela- 
tivistic outflow wi th a Lorentz factor of F ^ few (ZBSll, 
iBloom et al.l 1201 lb . The spectral energy distribution also 
demonstrated that the lack of detected optical variability re- 
quired s ignificant rest-frame extinction (Ay > 5 mag; ZBSl 1, 
iLevan et al. 2011), and that the X-rays were produced by 
a distinct emission component, rather than inverse Comp- 
ton scattering by the radio-emitting relativistic electrons 
(ZBSll). Finally, the evolution of the radio emission on a 
timescale of (5? 5 - 22 d pointed to an ambient density with 
a radial profile of roughly p oc r"^, as well as a mild increase 
in the energy of the outflow (ZBSl 1). 

The formation of a relativistic jet with dominant X-ray and 
radio emissi on were not predicted in standard tidal disruption 
models (e.g.. lReesll9"88tfstrubbe & Ouataertl20"09l) . which in- 
stead focused on the thermal optical/UV emission from the 
long-term accretion of the stellar debris. A signature of the 
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latter process is a mass accretion rate that evolves as M oc 
t~^l^, presumably lea ding to emission with the same tempo- 
ral dependence (e.g., Komossa & Greiner' 19 991: iGezari et alJ 
12008; van Velzen et al. 2011). Shortly before the discovery 
of Sw 1644+57. iGiannios & Metzgeil (1201 II) investigated the 
potential signature of a putative relativistic outflow, and con- 
cluded that the interaction of the outflow with the ambient 
medium will lead to radio emission on a timescale of 6t ^ \ 
yr (for typical off-axis observers). While the mechanism for 
the radio emission from Sw 1644+57 is interaction with an 
external medium, the actual light curves differ from the off- 
axis prediction. Toaddress this issue, in a follow-up pa- 
per iMstigiEiliD (1201 Ih (hereafter, MGMll) reconsidered 
the model for a relativistic jet interacting with an ambient 
medium. They draw on the inferences from the early radio 
emission described in ZBS 1 1 to infer the properties of the en- 
vironment and the jet kinetic energy, and use this information 
to predict the future evolution of the radio emission. 

This long-term radio evolution is of great interest because 
it can provide several critical insights: 

• The integrated energy release in the relativistic outflow, 
including the anticipated injection from on-going ac- 
cretion. 

• The density profile around a previously-dormant 
SMBH on 0.1 - 10 pc scales, which cannot be oth- 
erwise probed in AGN. 

• The potential to spatially resolve the outflow with very 
long baseline interferometry (VLBI), and hence to mea- 
sure the dynamical evolution (expansion and potentially 
spreading) of a relativistic jet. 

• Predictions for the radio emission from tidal disruption 
jets as viewed by off-axis observers on timescales of 
months to years to decades. 

The energy scale and jet dynamics are of particular impor- 
tance since the total energy input and the structure of the jet 
may also have implications for relativistic jets in GRBs and 
AGN. The ability to trace the environment on parsec scales 
provides a unique probe of gas inflow or outflow around an 
inactive SMBH on scales that cannot be probed outside of 
the Milky Way. Finally, the long-term radio emission from 
Sw 1644+57 will inform future radio searches for tidal dis- 
ruption events (TDEs) that can overcome the low detection 
rate in 7-rays/X-rays (due to beaming), and obscuration due 
to extinction in the optical/UV (as in the case of Sw 1644+57). 

To extract these critical properties we are undertaking long- 
term monitoring of the radio emission from Sw 1644+57 us- 
ing a wide range of centimeter- and millimeter-band facilities. 
Here we present radio observations of Sw 1644+57 that ex- 
tend to (5f « 216 d, and use these observations to determine 
the evolution of the total energy and ambient density. We find 
that the evolution of both quantities deviates from the behav- 
ior at (5f < 1 month (presented in ZBS 11), thereby providing 
crucial insight into the structure of the relativistic outflow and 
the ambient medium. This paper is the first in a series that 
will investigate the long-term radio evolution of Sw 1644+57 
and the implications for relativistic jets and parsec-scale envi- 
ronments around supermassive black holes, including efforts 
to resolve the source with VLBI and to measure polarization. 

The current paper is organized as follows. We describe the 
radio observations in ^ and summarize the radio evolution 



at 6t « 5-216 d in ^ In ^we present our modeling of the 
radio emission, which utilizes the formulation of MGMll. 
The implications for the energy scale and ambient density are 
discussed in 34.11 and 34.21 respectively, and we finally con- 
sider the implications for relativistic jets and the parsec-scale 
environments of SMBHs in ^ 

2. RADIO OBSERVATIONS 

Although Sw 1644+57 first triggered the Swift Burst Alert 
Telescope on 2011 March 28.55 UT, discernible 7-ray 
emission was detec ted starting on 2011 March 25 UT 
dBurrows et alj 120111) . We therefore consider 2011 March 
25.5 UT to be the actual initial time for the event and the 
associated relativistic outflow that powers the observed radio 
emission. Our radio observations of Sw 1644+57 commenced 
on 201 1 March 29.36 UT (19.4 hr after the Swift trigger and 
about 3.9 d after the initial 7-ray detection). Observations ex- 
tending to (5f « 26 d were presented in ZBS 1 1 . Here we report 
observations extending to w 216 d. Throughout the paper 
we use the standard cosmological constants with Hq = 70 km 
S-' Mpc"', Q.,„ = 0.27 and f^A = 0.73. 

2. 1 . Expanded Very Large Array 

We observed Sw 1644+57 with the EVLA** using the new 
Wideband Interferometric Digital Architecture (WIDAR) cor- 
relator to obtain up to 2 GHz of bandwidth at several frequen- 
cies. At all frequencies, we used 3C286 for bandpass and 
flux calibration. At 1.4 GHz, we used J 1634+6245 for phase 
calibration. For phase calibration at all other frequencies, 
we used J1638+5720, and also included a third calibrator, 
J1639+5357, at 5.8 GHz. The data were reduced and imaged 
with the Astronomical Image Processing System (AIPS) soft- 
ware package. The observations are summarized in Table [1] 
Minor changes to data values with respect to those reported 
in ZBS 1 1 are due to additional flagging of the data. The er- 
rors reported in Table [T] are statistical uncertainties only; the 
overall systematic uncertainty in the flux calibration is < 5%. 

2.2. AMI Large Array 

We observed with the AMI Large Array at 15.4 GHz with a 
bandwidth of 3.75 GHz. The maximum baseline is about 110 
m, with a resulting angular resolution of 25 arcsec. Observa- 
tions ranged in duration from 45 min to 1 1 hr. Observations of 
the compact source J1638+5720 were interleaved at intervals 
of 10 min as a phase reference, and the flux density scale was 
established by regular observations of the calibrators 3C48 
and 3C286. The telescope measures linearly-polarized signals 
(Stokes I+Q). The observations are summarized in Table[T] 

2.3. Combined Array for Research in Millimeter Astronomy 

We observed Sw 1644+57 with CARMA at frequencies of 
87.3 and 93.6 GHz with a total bandwidth ranging between 
6.8 and 7.8 GHz. We used Neptune as our primary flux cali- 
brator, and J1824+568 and J1638+573 as bandpass and phase 
calibrators. The overall uncertainty in the absolute flux cal- 
ibration is « 20%. Data calibration and imaging were done 
with the MIRIAD and AIPS software packages. The obser- 
vations are summarized in Table [T] We note some changes 

^ The EVLA is operated by the National Radio Astronomy Observatory, a 
facility of the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. The observations were obtained as part 
of programs lOC-145, llA-262, and llA-266 
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with respect to the flux densities reported in ZBSll, espe- 
cially those epochs where the total integration time on source 
was less than ^ 20 min (St « 10-25 d). 

2.4. Submillimeter Array 

We observed Sw 1644+57 with the SMA using at least 
seven of the eight antennas, in a wide range of weather condi- 
tions, with T225 ranging from 0.04 to 0.3. In each observation 
we combined the two sidebands, each with a bandwidth of 
4 GHz separated by 10 GHz, to increase the signal-to-noise 
ratio. The data were calibrated using the MIR software pack- 
age developed at Caltech and modified for the SMA, while for 
imaging and analysis we used MIRIAD. Gain calibration was 
performed using J1642+689, 3C345, and J1849+670. Abso- 
lute flux calibration was performed using real-time measure- 
ments of the system temperatures, with observations of Nep- 
tune to set the overall scale. Bandpass calibration was done 
using 3C454.3, J1924-292, and 3C279. The observations are 
summarized in Table [1] 

2.5. Very Long Baseline Interferometry 

We observed Sw 1644+57 with the NRAO Very Long Base- 
line Array (VLBA^) and the 100-m Effelsberg telescope in 
six epochs'" between 2011 April 2 and September 17 UT 
at 8.4 and 22 GHz. The observations were performed with 
eight frequency bands of 8 MHz bandwidth each in dual cir- 
cular polarization, r esulting in a to tal data rate of 512 Mbps. 
ICRFJ1638+5720 ( Ma et al.| [T998^, located only 0.92° from 
Sw 1644+57 was used for phase-referencing at both frequen- 
cies. At 22 GHz, we switched between the target and cal- 
ibrator every 40 s, while at 8.4 GHz we spent 40 s on the 
calibrator and 90 s on the target. This resulted in a total inte- 
gration time on Sw 1644+57 of about 93 min at 8.4 GHz and 
78 min at 22 GHz in epochs 1, 2, 5, and 6. Epochs 3 and 4 
were 1 hour shorter, resulting in integration times on source 
of 75 min at 8.4 GHz an d 60 min at 22 GHz. A second cali- 
brator, J 1657+5705 (Beas lev et 311 12002). was observed for 6 
min at each frequency to check the overall data calibration. At 
22 GHz, we also performed ~ 30-min blocks of geodetic ob- 
servations to perform atmospheric calib ration (for details see 
iBrunflialer et"ai]|2005HReid et al.ll2009h . 

The results from epoch 1 were presented in ZBSl 1. Here, 
we only discuss the 22 GHz data; the full VLBI data set 
will be presented together with future VLBI observations in 
an upcoming paper The data were correlated at the VLBA 
Array Operations Center in Socorro, New Mexico and cal - 
ibrated using AIPS and ParselTongue dKettenis et alj|2006h . 
We applied the latest values of the Earth orientation parame- 
ters and performed zenith delay corrections based on the re- 
sults of the geodetic block observations. Total electron con- 
tent maps of the ionosphere were used to correct for iono- 
spheric phase changes. Amplitude calibration used system 
temperature measurements and standard gain curves. We 
performed a "manual phase-calibration" using the data from 
one scan of J1638+5720 to remove instrumental phase offsets 
among the frequency bands. We then fringe-fitted the data 
from J1638+5720. Since J1638+5720 has extended struc- 
ture, we performed phase self-calibration, and later ampli- 
tude and phase self-calibration on J1638+5720 to construct 
robust models of J1638+5720 at both frequencies. These 

' The data were collected as part of programs BS210 and BS212. 

The dates of the observations are: 2011 April 2, April 8, May 5, May 
28, July 17, and September 17 UT. 



models were then used to fringe-fit the data again. Finally, the 
calibration was transferred to Sw 1644+57 and J1657+5705. 
The data were imaged in AIPS using robust weighting (with 

ROBUST=0). 

A linear fit to the positions of Sw 1644+57 from epochs 2-6 
gives a proper motion of 81 ± 37 /ias yr"' in right ascension 
and 46 ± 59 /zas yr"' in declination, with an additional sys- 
tematic uncertainty of about 40 /ias yr"'. This is consistent 
with no detectable motion. The corresponding 3(7 upper limit 
on the motion at the redshift of Sw 1644+57 is < 1 pc. Since 
the first epoch was correlated at a different position, it was not 
included in this fit to avoid second order systematic errors. 

To obtain an accurate estimate of the angular size in our lat- 
est epoch (2011 September 17 UT), we fit Gaussian models 
directly to the visibility data. We first self-calibrated the data 
in phase using solution intervals of 20 min, and then fitted 
a circular Gaussian model using a weighted least-squares fit. 
The source is not resolved in our observation. For unresolved 
sources the fitted size can be strongly correlated with the an- 
tenna amplitude gains, which are only imprecisely known. 
Our uncertainty on the fitted FWHM size, and consequently 
our upper limit, was estimated from a Monte-Carlo simulation 
in which we randomly varied the amplitude gains of the an- 
tennas by conservative factors of ±25% (for a more elaborate 
discussion of the uncertainties in estimating angular source 
sizes from a similar procedure see Bietenholz et al. 2010). 
The 3(7 upper limit on the FWHM source size is 0.22 mas, 
corresponding to a projected physical diameter of < 1.1 pc. 
For spherical expansion this corresponds to an upper limit on 
the apparent expansion velocity of < 3.8c. For a collimated 
relativistic source with an opening angle 9j the projected di- 
ameter is « AT^ctOj, indicating that for 6j = 0.1 the limit on 
the size of Sw 1644+57 coiTesponds to T < 5 at (5f ss 176 d. 
We compare this result with estimates of the source size from 
modeling of the radio emission in ^4.11 

3. OBSERVED EVOLUTION OF THE RADIO EMISSION 

The light curves at observed frequencies of 4.9, 6.7, 15.4, 
19, 24, 43, 90, and 230 GHz are presented in Figure □ As 
described in ZBSl 1, the light curves at 4.9 to 24 GHz exhibit 
an initial rapid increase, with a flux density of /v oc f'-^-f^, 
followed by a shallower increase, with oc f"^. In the mil- 
limeter band, the light curves peak on a timescale of a few 
days and then decline as /v oc f"'. This evolution is due to 
a synchrotron spectrum that is optically thin in the millime- 
ter band and optically thick in the centimeter band, with a 
peak frequency and flux density that decline as a function of 
time as Vp cx f"' -' and /v_p cx f""** (ZBSll). The properties 
of the spectral energy distribution, along with measurements 
of interstellar scintillation, established that the outflow is rel- 
ativistic (F K, 2-5), that the density profile of the circum- 
nuclear medium (CNM) is roughly pcNM r~^, and that the 
post-shock energy increases as £ oc /"'^ (ZBSl 1). 

The most striking result from our new observations is that 
the millimeter flux on a timescale of '-^ 100 d is significantly 
brighter than expected from an extrapolation of the early de- 
cline. Indeed, at 90 GHz the brightness is comparable to the 
initial peak aiSt ^ 10 d. A similar effect is observed at lower 
frequencies, where the light curves exhibit an upturn starting 
at « 30 d, with /v oc f"-^. The 90 GHz light curve is already 
in a declining (i.e., optically thin) phase at 5t > 100 d, while 
the Hght curves at 4.9 to 24 GHz peak at 5f « 100 d (19 and 
24 GHz), St « 125 d (15.4 GHz), St « 170 d (8.4 GHz), and 
(5f « 190 d (4.9 and 6.7 GHz), with peak flux densities of w 3 1 
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mJy, « 28 mJy, « 23 mJy, and « 20 mJy, respectively. These 
observations are at odds with the light curve evolution pre- 
dicted by MGMl 1 based on our radio data at 5 -22 d, using 
a model that assumes a constant energy and a steady p oc 
CNM profile (Figure [B- 

The observed change in the light curve evolution at 5t > 
1 month requires an increase in the outflow kinetic energy 
and/or the CNM density. The evolution of the peak frequency, 
roughly i/p oc f"' at 6t ~ 30- 100 d, is slower than expected for 
the characteristic synchrotron frequency, iy,„ oc (see ®. 
Since Vm only depends on the kinetic energy, v,,, oc Zs^^f"^/^ 
(iGranot & S ari 2002), the observed evolution requires an in- 
crease in the energy, with Ej (xt. At the same time, the similar 
light curve evolution at 4.9 and 24 GHz, which straddle the 
self-absorption frequency (see ®, indicates that the density 
must also change. In particular, at 4.9 GHz we expect /v cx 
£j.5^-i.5 ^jjjjg 24 GHz we expect oc Ef^A°-^ f-*'-^; 
here A is a fiducial density that parametrizes the CNM den- 
sity, with p = Ar~^. For the observed light curve evolution of 
/v oc f "-^ at both frequencies, and using the inferred evolution 
of Ej, we find that A oc t^-^^. Thus, both the energy and den- 
sity normalization appear to increase at St > 30 d. In the next 
section we model the radio emission in detail to extract the 
temporal evolution of the energy and density. 

4. DETAILED MODELING OF THE RADIO EMISSION 

To determine the temporal evolution of the relativistic out- 
flow and the radial density profile we use the formulation 
of MGMl I, w hich draws on the GRB afterglow model of 
IGranot & Sari ( 12002.) (hereafter, GS02). The MGMl 1 model 
is largely motivated by the basic properties of the outflow 
that were determined in our initial study of the radio emis- 
sion (ZBSll), as well as b y the observed evolution of the X- 
ray emission (Burro ws et al.ll201 ID . The model assumes that 
the outflow is collimated, with an opening angle 9j, and has a 
Lorentz f actor of F , ; in ZB S 1 1 we demonstrated that Tj ~ few 
(see also ' Bloom et al.ir201 Ih . The kinetic luminosity of the 
outflow, Lj.iso, is assumed to be constant for a timescale tj, fol- 
lowed by a decline as Lj^so oc f"^''-' at f > tj. This is expected 
for fallback accretion in a TDE, and is observationally mo- 
tivated by the evolution of the X-ray light cu rve, which also 



indicates that f , w 10^ s ( Burrows et al.ll20rT i). As a result of 



the long-term injection of energy, the expectation is that the 
kinetic energy will gradually increase from an initial value of 
£i,iso = Lj.isotj to a final level of (5/2)Lj isof). From the initial 
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« 6 X 10 erg s , the expected total 
3 X 10^-'(ex/0.5)"' erg, where ex is the ef- 



X-ray luminosity, Lx> 
energy is £j j 

ficiency of the jet in producing X-rays (MGMll), while the 
predicted beaming-corrected energy is isy « 1 .5 x 10^' erg. 

The interaction of the outflow with the CNM leads to syn- 
chrotron emission due to the acceleration of electrons and am- 
plification of magnetic fields. Based on the early radio obser- 
vations, ZBS 1 1 demonstrated that the CNM density roughly 
follows pcNM oc (hereafter. Wind medium), which is the 
profile adopted for the analysis in MGMll. The Lorentz fac- 
tor of the fluid behind the forward shock is Fsh.ps = - 
T j^lncHMn nj), where mcnm is the circumnuclear density 
and rij = Lj ifio/ATTr^mpC^T^j is the density of the ejecta (see 
MGMl 1); this relation is appropriate for the case of a Newto - 
nian reverse shock, i.e., Fy <C «j/ncNM (Sari & Piran 1995!), 
which as we find below is appropriate for Sw 1644H-57 (Ta- 
ble |2|i. A key difference between the MGMl 1 model and the 



GS02 GRB afterglow model is that the former includes a sup- 
pression of the flux density by a factor of (Fsh.Fs6'j)^/2 due to 
the finite extent of the collimated outflow (i.e., Fsh,FS ^ 1 / ^j)- 
Thus, the synchrotron emission from the relativistic outflow 
is determined by Lj iso, hcnm, and the fractions of post-shock 
energy in the radiating electrons (e^) and the magnetic field 

We use this model (Equations 7-1 1 of MGMl 1), along with 
the smoothing formulation of GS02, to fit snap-shot broad- 
band spectra of Sw 1644H-57 on a timescale of St w 5-216 
d. In each epoch we fix = eg = 0.1 and p = 2.5, and we 
then determine the best-fit values of Lj iso and hcnm; here p 
is the power law index of the electron Lorentz factor distribu- 
tion, A^(7) oc 7"''. For our choice of p value the synchrotron 
frequencies (self-absorption: I'a', peak: t/,,,) and flux normal- 
ization [Fi^iv = Va)] for t > tj are given by (MGMl 1, GS02): 



iyJt) = 42xl(fe^l,e^l,Lr-^,,tj},n 
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z.„(f) = 3.5 X lQ''e^^AU^ioA,t]l [-} Hz, (2) 
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^;[i^(f)] = 345e,,_iLiiL,48f^6«Tr^: 
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where we use the notation X = WXy, and nig is the CNM 
density (hcnm) at a fiducial radius of 10'** cm. At f < tj we 
use the same equations but with the time dependencies mod- 
ified to OC t~\ u,n OC f"', and F^iv = Ua) oc (MGMll). 
The synchrotron sp ectrum is given by (c.f.. Spectrum 1 of 
IGranot & Sarill2002h : 



1+1^ 



-Si A 



-lAi 



l2(/32-A)' 



(4) 



where si and S2 are smoothing parameters (GS02), and /3i = 2, 
= 1/3, and (3^ = (l-p)/2 are the power law indices for each 

segment of the synchrotron spectrum. 
Equations [T] and |4] are appropriate when i^a < Vm, and have 

to be modified when v„i < Va as foflows (GS02, MGMl 1 "): 



Va OC / 



(5) 



fix 






+ 



















.!2(/32-A)' 



-1/.S2 



(6) 



where the spectrum is now normalized at v = iy„„ s-i and S2, 
take on different values (GS02) and /3i = 2, (32 = 5 /2, and 
A = (1 -p)/2. To smoothly connect the evolution in the two 
phases we use a weighted average of Equations |4] and |6] with 
the weighting determined by the time difference relative to 
the transition time, 4^, defined as i'a(ram) = i^mitam)- We find 

' ' Note that MGMl 1 use the scaHngs appropriate for Spectrum 3 of GS02, 
while we use the appropriate Spectrum 2 with i/,,, < !/„ < Uc. 
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that fam « 275 d, so this transition does not affect the data 
presented in this paper 

Instead of imposing a specific temporal evolution on Lj iso 
and ncNM, we model each broad-band radio spectrum inde- 
pendently to extract the evolution of these quantities, and in 
turn the time evolution of the emission radius (r) and Fsh (and 
hence Tj), as well as the radial density profile. The results 
of the individual fits are shown in Figure |2] and the relevant 
extracted parameters are plotted in Figure [3] and summarized 
in Table|2] The broad-band SEDs reveal a complex evolution. 
At (5f « 5 - 22 d both the peak frequency and peak flux den- 
sity decrease with time (as noted by ZBSll), but starting at 
(5f « 36 d the peak flux density begins to increase, reaching a 
maximum at St ^ 100 d, and subsequently declining again. 

The VLBI size Hmit of rproj < 0.55 pc at St « 175 d 
can be used to set an independent upper bound on the ratio 
Lj.isoA»/n\s. In particular, on this timescale we expect rproj ~ 
2.9t j^(,(Lj,i^oMi/nis.)°-^9j pc, indicating that Lj.iso,48/ni8 < 3.6 
(for 6j = 0.1). From the results of our radio modeling we find 
on a similar timescale that -Lj,iso,48 ~ 0.45 and «ig « 3 (Ta- 
blelU, and hence Lj_iso,48/ni8 ~ 0.15 is in agreement with the 
VLBI size limit. On the other hand, if we assume = 0.01 
the resulting value of Lj_iso,48 increases by about a factor of 
15, while «i8 decreases by about a factor of 2, leading to 
io,iso,48/ni8 ~ 4, in violation of the VLBI size limit'^ We 
therefore conclude from the VLBI results that eg > 0.01. 

4. 1 . The Evolution ofEj iso 

A detailed analysis of the evolution of i',,, (shown in Fig- 
ure |4]i reveals a complex behavior. While it continuously de- 
clines as a function of time, the decline rate is always shal- 
lower than the expected f"^/^ for constant energy. Since i',,, 
only depends on Lj i^o and the equipartition fractions (which 
are not expected to change with time), the shallower than ex- 
pected evolution directly implies that Lj iso continuously in- 
creases with time. In particular, at 5-22 d and > 100 d, the 
observed evolution of ly^ oc f"'-'^ indicates that Lj iso oc f"'^, 
while in the intermediate phase (22-100) d, i^,,, oc f"''-^^ points 
to a more rapid increase in the energy scale, with Lj iso oc f' ' . 

The time evolution of Lj iso is shown in Figure |5] We 
find that at St < 20 d Lj.iso « 4.5 x lO''^ erg s"' (or fij.iso « 
4.5 X 10''^ erg for tj = 10^ s). Changing the values of the 
equipartition fractions from our assumed values = = 0.1, 
we find that Lj iso increases by about a factor of 3 if = 0.01 
or by a factor of 15 if = 0.01 {eg < 0.01 is ruled out by the 
VLBI limits on the source size). These scalings hold for the 
overall time evolution of Ljiso- 

Also shown in Figure q\ is the expected evolution for a 
model with Lj iso = const at f < tj and Lj ^so cx: t~^^^ at t > tj, 
as inferred from the X-ray light curve (iBurrows et al.ll20Tll 
MGMll). In this model, we expect the integrated value of 
Lj iso (i-C-, the total isotropic-equivalent energy) to approach 
(5/2)Lj,iso as f — > oo, and to roughly double within St ~ 5tj ~ 
50-60 d. This seems to be the case based on the radio data 
at St < 30 d, but the subsequent rapid increase in energy by 
about an order of magnitude at 30-216 d cannot be explained 
with continued injection from a f"^/^ tail. On the other hand, 
the X-ray light curve of Sw 1644H-57 (Figure |5]l agrees well 
with this simple luminosity evolution, with Lx,iso,o ~ 6 x 10'*^ 
erg s"' at f < 13 d. A comparison to the inferred value of 

Changing to eg = 0.01 has almost no effect since both the energy and 
density increase by about a factor of 3 compared to eg = 0.1. 



Lj indicates a high efficiency''' of producing X-rays of 
ex « 0.5. 

Since continuous energy injection from an L oc f tail 
cannot explain the inferred rise in integrated kinetic energy, 
the radio data require a different energy injection mechanism. 
One possibility is that the outflow has a distribution of Lorentz 
factors, with increasing energy as a function of decrea sing 
Lorentz factor, Ej(> Tj) oc Pj"" (ISari & Meszarosll2000l ). In 
this scenario the energy will increase with time as material 
with a lower Lorentz factor catches up with decelerated ejecta 
that initially had a higher Lorentz factor. The continuous in- 
jection of energy will also lead to a more rapid increase in 
radius, with r oc f(i+^)/(7+.'-2"i) (Sari & Meszaros||2000|); here 
m is the power law index describing the CNM density profile 
(p oc r"™). In Figure[3]we plot the inferred radius as a function 
of time, and find that it follows r oc f"^, faster than expected 
in a simple Wind model (r oc f°-^). We also find from the den- 
sity profile that m w 1 .5 (see 34.2l i. thereby leading to w 3.5. 
Thus, we expect in this scenario a fairly steep profile for the 
energy as a function of Lorentz factor of E{> Tj) oc Fy^^. 

With this relation, and the inferred time evolution of the jet 
Lorentz factor (FigureO we can infer the expected increase in 
energy. We find F^ « 5.5 at St = 10 d and F^ « 2.2 at St = 216 
d, indicating an expected increase in energy by about a factor 
of 10 for our inferred profile. This is in excellent agreement 
with the data. Thus, a jet with a distribution of Lorentz fac- 
tors of E(> Tj) oc Fy^^ naturally explains the evolution of 
source size and the substantial increase in energy beyond the 
on-going input from a t~^^^ tail. 

4.2. The Radial Density Profile 

In addition to the substantial increase in energy, we also find 
that the normalization of the density profile, nig = ncNM'"?8' 
changes as a function of time (Figure O. This indicates that 
the radial density structure is not a simple Wind profile, for 
which ni8 is by definition constant. In particular, we find that 
«i8 decreases at 5-10 d, followed by a steady increase at 
10-216dis about a factor of 8. The resulting radial density 
profile is shown in Figure |6] and clearly reflects the complex 
evolution of The density profile at r w 0.1 -0.4 pc is p oc 
r"'^, followed by a uniform density at r « 0.4-0.6 pc, and a 
slow transition back to p oc by r w 1.2 pc. The inferred 
density is « « 0.2-2 cm"^ at r w 1.2-0.1 pc, respectively. 
These values scale with the choice of equipartition fractions: 
for es = 0.01 the density is about 3 times larger, while for 
ee = 0.01 it is about 1.7 times lower Thus, the inferred density 
is fairly robust to large changes in the equipartition fractions. 

The flattening at r « 0.4-0.6 pc coincides with the rapid 
increase in energy at (5f ~ 35 - 100 d. This can be understood 
in the context of an outflow with a distribution of Lorentz 
factors since the relative increase in density compared to the 
previous r~'^ profile leads to enhanced contribution from 
lower Lorentz factor ejecta. A similar effect has been ob- 
served in the radio emission from some core-collapse super- 
novae, which have a steep ejecta profile with E oc (e.g. , 
ChevalieiifT98llMatzner & McKeelfT99l iBerger et aLllBol 
Soderberg et aLll2006l) . 

Also shown in Figure|6]are the inferred densities in the cen- 
tral parsec of the Milky way galactic center from X-ray mea- 
surements ( Baganoff et al. 2003). The Galactic Center den- 
sity is similar that inferred for Sw 1644H-57 at r w 0.05 pc, but 

The X-ray efficiency is instead ex ~ 0.3 if eg = 0.01 or 0. 1 if = 

0.01. 
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is about 30 times larger at r w 0.4 pc. If the bulk of the gas 
in th e central parsec is due to mass loss from massive stars 
(e.g.. lMelia|[T99l iBaganoff et al.1 120031; IOuataerlll2004 ). the 
lower density inferred here may indicate a much smaller num- 
ber of massive stars in the central parsec of the host galaxy of 
Sw 1644+57. Although we cannot investigate star formation 
on scales smaller than ~ 1 kpc in the host of Sw 1644+57, 
the overall star formation rate in this galaxy is indee d a fac- 
tor of several times lower than in the Milky Way ( Leva n et aP 
1201 iV Regardless of the reason for the difference in density, it 
is quite remarkable that the radio emission from Sw 1644+57 
offers as detailed a view (or better) of the density profile in the 
inner parsec around an inactive SMBH at z = 0.354 as avail- 
able for the Galactic Center at a distance of only 8.5 kpc. 

4.3. Optical and Near-Infrared Emission 

Using the radio modeling we can also predict the optical 
and near-IR emission from the jet. In Fig ure|7lwe plot the op- 
tical r-band upper limits froml Levan et al. (2011), as well as 
their near-IR A'-band measurements. The latter include the to- 
tal flux from Sw 1644+57 and its host galaxy since transient- 
free templates are not presently available. We estimate the 
host contribution to be about 20 /iJy (K « 20.6 AB mag) 
by requiring that the overall shape of the /T-band light curve 
match our predicted light curve. We note that this only af- 
fects the light curve shape at St > 15 d since at earlier times 
the observed flux is dominated by the transient itself. A sim- 
ple extrapolation of our model over-predicts the /iT-band flux 
density by about an order of magnitude, indicating the pres- 
ence of an additional break in the spectrum between the ra- 
dio and optical/near-IR band. Such a break is indicative of 
the synchrotron cooling frequency, Vc- To explain the /T-band 
flux density requires i^c ^ 10'^ Hz. Even if we include this 
break, the observed r-band limits are still a factor of about 70 
times fainter than the model prediction. To explain this dis- 
crepancy reguiresjiostgalaxy extinction of Ay.host ^ 3.5 mag 
(see also iLevanetaTI 120 111 ZBSll). With the addition of a 
cooling break and host galaxy extinction we find an excellent 
match between our model and the near-IR fluxes. Since these 
data were not used in the fitting, the agreement between our 
model and the data provides an independent confirmation of 
the model. 

5. IMPLICATIONS FOR RELATIVISTIC JETS AND THE 
ENVIRONMENTS OF SUPERMASSIVE BLACK HOLES 

Our on-going radio observations have uncovered two unex- 
pected and critical results regarding the nature of the relativis- 
tic outflow from Sw 1644+57 and the parsec-scale environ- 
ment around a previously-dormant SMBH. We find that the 
relativistic outflow is not dominated by a single Lorentz fac- 
tor, as typically assumed in GRB jet models. Instead, the out- 
flow has a Lorentz factor distribution with Ej{> Tj) oc 
at least over the range of Lorentz factors probed so far, Tj w 
2.2-5.5. This scenario is reminiscent of supernova ejecta, 
in whi ch the energy is a strong function of velocity, E cx v"^ 
dChev alier 1982; Matzner & McKee 1999). In the SN case, 
this coupling reflects shock acceleration through the steep 
density gradient in the outer envelope of the star In the case 
of Sw 1644+57 we do not expect such a density gradient, sug- 
gesting instead that the ejecta structure may be an intrinsic 
property of relativistic jet laun ching, potentially thro ugh the 
Blandford-Znajek mechanism (Bla ndford & Znaje k 1977). In 
this case, the same structure may also apply to the relativistic 
jets of GRBs, which are generally assumed to have a single 



Lorentz factor. 

To assess this possibility we compare the ratio of 
event duration to the accretion disk dynamical timescale, 
fdyn oc M-Rnir/ rgf'l^, where rg = IGM^u/c^ is the black 
hole Schwarzschild radius. For similar ratios of (r/rg), 
since Mbh,Sw1644+57 ~ IO^'Mq - IO^Mbh.grbs, the dynami- 
cal timescale in Sw 1644+57 is about 10^ times longer than 
in GRBs. The ratio of durations is similar, ^ 10'' s for 
Sw 1644+57 and - 10 s for GRBs. Therefore, the ratio of 
duration to dynamical timescale is similar, and the jet launch- 
ing mechanism may imprint a similar profile in the case of 
Sw 1644+57 and GRBss. Of course, in long GRBs the jet 
still has to propagate through the stellar envelope, which may 
modify the jet structure. 

We next address the potential implications of the radial den- 
sity profile and the flattening at r w 0.4-0.6 pc. One possi- 
bility is that this scale represents the Bondi radius for spheri- 
cal accretion from the CNM, rs = GM/c^ = (c/cs)^ here 
Cj is the sound speed in the CNM and rg = 2GM/c^ is the 
Schwarzschild radiu s of the SMBH. For a CNM gas tem- 
peratu re of - 10^ K (Bag anoff et al.ll2003l: lNaravan & FabianI 
1201 lb . we find rg « lO^r, « 0.1 - 1 pc for the'SMBH mass 
of lO''- 10^ M0 inferred for Sw 1644+57. Thus, the radius 
at which we infer a density enhancement is in reasonable 
agreement with the expected Bondi radius. In addition, we 
note that the expected density profile inside the Bondi radius, 
p cx r"^/^, matches the inferred profile for Sw 1644+57 on a 
scale of ^ 0.15-0.4 pc. In the standard scenario, the medium 
outside the Bondi radius is assumed to have a uniform den- 
sity, while here we infer a continued decline in the density at 
r w 0.6- 1.2 pc. This may reflect the complex conditions in 
the inner few parsecs around the SMBH in which the inter- 
stellar density profile may be mainly influenced by mass loss 
fro m stars (as in the ca se of the Milky Way galactic center; 
e.g.. lKrabbe et alJll991h . Indeed, from the inferred density of 
"CNM ~ 0.6 cm"-^ at r w 0.5 pc, the required mass loss rate is 
M w 5 X 10"^ Mq yr"' for a wind velocity of lO-' km s"' . 

In particular, in the case of the Galactic Center it has long 
been believed that the gas in the central parsec is supplied 
by mass loss from massive stars, with the bulk of the gas 
(> 90%) being thermally expelled in a wind (Melia 19921 
Baganoff et al. 2003; Ouataert 2004, but see Sazonov et al] 
20 111 for an alternative interpretation of the X-ray emission). 
Using a spherically symmetric hyd rodynam i c simu lation of 
the gas supplied by stellar winds, lOuataerlj (l2004l) showed 
that beyond about 0.4 pc gas is mainly expelled with a result- 
ing Wind profile (p cx r ) on larger scales, while on scales 
of < 0.2 pc the profile is roughly as expected for Bondi ac- 
cretion (p cx r"-'/^); see Figure |6l The profile we find here is 
somewhat different from this model, but this may be due to 
a different distribution of massive stars, or to the simplifying 
assumption of spherical symmetry in the model. 

6. CONCLUSIONS AND FUTURE DIRECTIONS 

We presented radio observations of Sw 1644+57 extending 
to St 216 d and spanning a wide range of frequencies. The 
evolution of the radio emission changes dramatically atSt>l 
month, requiring an increase in the total energy by about an 
order of magnitude, a density profile of p cx r"^/^ (0.1 - 1.2 
pc), and a flattening at r « 0.4-0.6 pc. A comparison of 
the model to optical limits and near-IR detections indicates a 
cooling break at i^^ ~ 10'^ Hz and host galaxy extinction of 
Av,ho,st > 3.5 mag. 
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The increase in energy cannot be explained by injection 
from an L c>c t~^l^ tail that is expected in tidal disruption events 
and which matches the evolution of the X-ray emission. We 
conclude that a natural explanation is a structured outflow 
with £■(> Tj) oc rj^^. The inferred density profile and the ra- 
dial scale of the density enhancement are in rough agreement 
with the expectation for Bondi accretion from a circumnuclear 
medium. The jet energetics and structure, as well as the de- 
tailed density profile on ^ 0.1 - 1 pc scale are a testament to 
the important insight that can be gained from continued radio 
observations of Sw 1644H-57. In particular, the radial density 
profile is traced in greater detail than even the inner parsec of 
the Milky Way. Continued radio observations will probe the 
environment to a scale of ^ 10 pc in the coming decade. 

Using the results of our analysis we can predict the fu- 
ture evolution of the radio emission (modulo any future un- 
predictable changes in energy and/or density as we have 
found here). We use the evolution of Lj iso and «cnm as 
inferred from the data at 5f < 216 d, and assume that the 
density will continue to evolve as p oc r"'-^ and that the en- 
ergy will increase to a maximum beaming-corrected value 
of =Lj,isoO[l-cos(6'j)] with Ej = \Q^^, 3 x 10^^ erg. 
The resulting light curves at 6 and 22 GHz are shown in Fig- 
ure [8] The long-term evolution is marked by a break when 
Ej achieves its maximum value, corresponding to about 5, 28, 
and 180 yr for our three choices of maximum energy. Us- 
ing the 5a sensitivity of the EVLA in an observation of a few 
hours'"*, we find that the emission at 22 and 6 GHz should be 
detectable for at least ^ 40 yr and ^ 80 yr, respectively. In- 
deed, any significant upgrades to the EVLA or the construc- 
tion of more sensitive radio facilities in the coming decades 
may extend the range of detectability to centuries'^. The same 
is true if the total energy scale is ^ 10^^ erg. 

An equally important question is whether the jet will be 
resolvable with VLBI in the future. The projected radius is 
Voj ~ '"^i' l^'^S J^*- maintains its collimation. In 

Figure |9] we plot the predicted future evolution of r using the 
prescription described above. We find that for 6j = 0.1 and 
a best-case VLBI angular resolution of « 0.2 mas (FWHM), 
the source should become resolvable at (5f « 6 yr On this 
timescale the 22 GHz flux density is expected to be only w 
2 mJy (Figure |8]l, still accessible with VLBI. While the flux 
density at 6 GHz is expected to be larger by about a factor 
of 2.6, the angular resolution at this frequency is poorer by 

At 22 GHz we use the sensitivity for the full 8 GHz bandwidth that will 
become available some time in 2012. 

Significant budget cuts to radio facilities in the future may lead to the 



about a factor of 3.7, making it less competitive than 22 GHz. 
Thus, we conclude that the radio emission from Sw 1644-1-57 
may be marginally resolved in a few years. On the other hand, 
if the jet undergoes significant spreading on the timescale at 
which it becomes non-relativistic (as expected for GRB jets: 
e.g., Li vio & Waxman 2000.) it is possible that it will become 
resolvable at 5/ ~ 1 - 2 yr when the expected 22 GHz flux 
density is still ~ 10 mJy. 

We are undertaking continued multi-frequency radio mon- 
itoring of Sw 1644-1-57 to follow the long-term evolution of 
the relativistic outflow and the radial profile of the ambient 
medium. Even in the absence of any future dramatic changes 
relative to the current evolution, we expect that in the next few 
years we may be able to determine the total energy of the rel- 
ativistic outflow, measure the spreading of the jet, and study 
the radial density profile to a scale of ~ 10 pc. Future papers 
in this series will detail these results. 
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TABLE 1 

Radio Observations of Sw 1644+57 



St" 


Facility 


Frequency 


Flux Density 


(d) 




(GHz) 


(mjy) 


6.79 


EVLA 


1.4 


0.21 ±0.08 


126.59 


EVLA 


1.4 


l.lOitO.lO 


174.47 


EVLA 


1.4 


1.60±0.11 


197.41 


EVLA 


1.4 


1.80±0.10 


3.87 


EVLA 


4.9 


0.25 ± 0.02 


4.76 


EVLA 


4.9 


0.34 ±0.02 


5.00 


EVLA 


4.9 


0.34 ± 0.02 


5.79 


EVLA 


4.9 


0.61 ±0.02 


6.78 


EVLA 


4.9 


0.82 ±0.02 


7.77 


EVLA 


4.9 


1.48 ±0.02 


9.79 


EVLA 


4.9 


1.47 ±0.02 


14.98 


EVLA 


4.9 


1.80 ±0.03 


22.78 


EVLA 


4.9 


2.10±0.01 


35.86 


EVLA 


4.9 


4.62 ±0.02 


50.65 


EVLA 


4.9 


4.84 ±0.03 


67.61 


EVLA 


4.9 


5.86 ±0.03 


94.64 


EVLA 


4.9 


9.06 ±0.03 


111.62 


EVLA 


4.9 


9.10±0.03 


126.51 


EVLA 


4.9 


9.10±0.03 


143.62 


EVLA 


4.9 


11.71 ±0.03 


164.38 


EVLA 


4.9 


12.93 ±0.05 


174.47 


EVLA 


4.9 


12.83 ±0.06 


197.41 


EVLA 


4.9 


13.29 ±0.03 


213.32 


EVLA 


4.9 


12.43 ±0.04 


3.87 


EVLA 


6.7 


0.38 ±0.02 


4.76 


EVLA 


6.7 


0.63 ± 0.02 


5.00 


EVLA 


6.7 


0.64 ±0.02 


5.79 


EVLA 


6.7 


1.16±0.02 


6.79 


EVLA 


6.7 


1.47 ±0.02 


7.77 


EVLA 


6.7 


1.50 ±0.02 


9.79 


EVLA 


6.7 


2.15 ±0.02 


14.98 


EVLA 


6.7 


3.79 ±0.03 


22.78 


EVLA 


6.7 


3.44 ±0.01 


35.86 


EVLA 


6.7 


6.39 ±0.02 


50.65 


EVLA 


6.7 


5.70 ±0.02 


67.61 


EVLA 


6.7 


8.94 ±0.03 


94.64 


EVLA 


6.7 


13.43 ±0.03 


111.62 


EVLA 


6.7 


13.66 ±0.03 


126.51 


EVLA 


6.7 


14.16 ±0.04 


143.62 


EVLA 


6.7 


16.85 ±0.04 


164.38 


EVLA 


6.7 


18.27 ±0.06 


174.47 


EVLA 


6.7 


19.59±0.16 


1 Q7 A 1 


FVT A 

Hy V ijrt 


D. / 




Zi J. JZ 


FVT A 
C V ijrt 


D. / 




14.97 


EVLA 


8.4 


5.49 ±0.09 


127.69 


EVLA 


8.4 


19.03 ±0.14 


159.77 


EVLA 


8.4 


22.15 ±0.20 


174.47 


EVLA 


8.4 


23.19±0.38 


177.50 


EVLA 


8.4 


23.65 ±0.16 


197.41 


EVLA 


8.4 


22.42 ±0.10 


213.32 


EVLA 


8.4 


22.04 ±0.13 


219.22 


EVLA 


8.4 


21.52 ±0.09 


5.81 


AMI-LA 


15.4 


2.69 ±0.44 


6.64 


AMI-LA 


15.4 


3.62 ±0.22 


7.62 


AMI-LA 


15.4 


4.32 ±0.28 


8.55 


AMI-LA 


15.4 


5.07 ±0.36 


9.56 


AMI-LA 


15.4 


6.68 ±0.51 


10.78 


AMI-LA 


15.4 


6.74 ±0.48 


11.56 


AMI-LA 


15.4 


7.50 ±0.44 


13.64 


AMI-LA 


15.4 


8.02 ±0.71 


14.57 


AMI-LA 


15.4 


8.43 ±0.32 


16.47 


AMI-LA 


15.4 


8.86 ±0.58 


18.65 


AMI-LA 


15.4 


8.62 ±0.40 


19.73 


AMI-LA 


15.4 


10.04 ±0.63 


21.78 


AMI-LA 


15.4 


10.91 ±0.89 


22.76 


AMI-LA 


15.4 


11.01 ±0.77 


25.38 


AMI-LA 


15.4 


10.28 ±0.58 


26.74 


AMI-LA 


15.4 


11.36 ±0.80 


31.45 


AMI-LA 


15.4 


11. 24 ±0.41 


33.61 


AMI-LA 


15.4 


12.14±0.16 


34.79 


AMI-LA 


15.4 


11. 89 ±0.25 


35.71 


AMI-LA 


15.4 


13.39 ±0.49 


37.49 


AMI-LA 


15.4 


13.63 ±0.45 


38.60 


AMI-LA 


15.4 


13.72 ±0.39 
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TABLE 1 — Continued 



St" 


Facility 


Frequency 


nux Density 






(GHz) 




39.78 


AMI-LA 


15.4 


10.80 ±0.91 


40.67 


AMI-LA 


15.4 


13.38 ±0.40 


41.69 


AMI-LA 


15.4 


13.64±0.21 


43.63 


AMI-LA 


15.4 


14.21 ±1.24 


45.74 


AMI-LA 


15.4 


12.33 ±0.37 


48.69 


AMI-LA 


15.4 


14.06 ±0.41 


49.64 


AMI-LA 


15.4 


13.94 ±0.05 


50.54 


AMI-LA 


15.4 


14.39 ±0.25 


53.39 


AMI-LA 


15.4 


15.94 ±0.59 


55.62 


AMI-LA 


15.4 


15.28 ±0.33 


56.52 


AMI-LA 


15.4 


17.91 ±0.07 


60.73 


AMI-LA 


15.4 


15.00±1.18 


62.54 


AMI-LA 


15.4 


19.23 ±0.62 


63.45 


AMI-LA 


15.4 


16.47 ±0.21 


65.51 


AMI-LA 


15.4 


18.77 ±0.37 


68.60 


AMI-LA 


15.4 


19.06 ±0.48 


70.69 


AMI-LA 


15.4 


20.34 ± 0.38 


71.38 


AMI-LA 


15.4 


19.36 ±0.56 


73.48 


AMI-LA 


15.4 


21.27 ±0.34 


74.48 


AMI-LA 


15.4 


21.84 ±0.59 


76.65 


AMI-LA 


15.4 


2 1.40 ±0.51 


77.47 


AMI-LA 


15.4 


23.08 ±0.26 


78.49 


AMI-LA 


15.4 


23.26 ±0.65 


79.49 


AMI-LA 


15.4 


23.22 ±0.24 


80.47 


AMI-LA 


15.4 


22.94 ±0.50 


81.48 


AMI-LA 


15.4 


21 .74 ±0.45 


83.42 


AMI-LA 


15.4 


23.99 ±0.50 


86.65 


AMI-LA 


15.4 


22.94 ±0.33 


88.37 


AMI-LA 


15.4 


24.69 ±0.26 


89.64 


AMI-LA 


15.4 


25.90 ±0.50 


91.58 


AMI-LA 


15.4 


26.07 ±1.04 


92.58 


AMI-LA 


15.4 


25.40 ±0.55 


95.39 


AMI-LA 


15.4 


25.34 ±0.72 


98.42 


AMI-LA 


15.4 


26.15 ±0.25 


100.61 


AMI-LA 


15.4 


27.83 ±0.57 


101.60 


AMI-LA 


15.4 


25.96 ±0.91 


102.60 


AMI-LA 


15.4 


26.60 ± 0.26 


105.59 


AMI-LA 


15.4 


26.68 ±0.21 


107.38 


AMI-LA 


15.4 


27.89 ±0.33 


108.32 


AMI-LA 


15.4 


26.84 ± 0.40 


110.34 


AMI-LA 


15.4 


28.24 ±0.22 


112.50 


AMI-LA 


15.4 


27. 10 ±0.58 


113.50 


AMI-LA 


15.4 


28.82±0.31 


115.49 


AMI-LA 


15.4 


28.39 ±0.23 


118.31 


AMI-LA 


15.4 


28.58 ±0.47 


119.48 


AMI-LA 


15.4 


26.90 ±0.80 


120.55 


AMI-LA 


15.4 


27.88 ±0.35 


121.55 


AMI-LA 


15.4 


26.92 ±0.73 


122.18 


AMI-LA 


15.4 


29.87 ±0.61 


123.45 


AMI-LA 


15.4 


28.86 ±0.41 


124.38 


AMI-LA 


15.4 


27.46 ± 0.97 


125.51 


AMI-LA 


15.4 


27.40 ±0.56 


127.38 


AMI-LA 


15.4 


28.96 ±0.39 


128.35 


AMI-LA 


15.4 


28.69 ±1.00 


130.38 


AMI-LA 


15.4 


27.87 ±0.64 


131.41 


AMI-LA 


15.4 


28.94±1.11 


132.49 


AMI-LA 


15.4 


29.39 ±0.96 


133.35 


AMI-LA 


15.4 


30.81 ±1.03 


134.36 


AMI-LA 


15.4 


29.73 ±0.52 


135.32 


AMI-LA 


15.4 


31.25 ±0.59 


136.46 


AMI-LA 


15.4 


29.31 ±1.29 


137.51 


AMI-LA 


15.4 


29.58 ±0.23 


138.51 


AMI-LA 


15.4 


28.50 ±0.69 


139.51 


AMI-LA 


15.4 


28.96 ±0.57 


140.50 


AMI-LA 


15.4 


29.22 ±0.52 


141.13 


AMI-LA 


15.4 


29.03 ± 0.45 


142.30 


AMI-LA 


15.4 


26.60 ±0.31 


143.26 


AMI-LA 


15.4 


28.96 ±0.61 


144.42 


AMI-LA 


15.4 


28.36 ±0.57 


146.48 


AMI-LA 


15.4 


28.92 ±0.64 


147.48 


AMI-LA 


15.4 


27.44 ± 0.76 


148.27 


AMI-LA 


15.4 


29.49 ±0.97 


149.27 


AMI-LA 


15.4 


29.82 ±0.43 


150.47 


AMI-LA 


15.4 


29.83 ±0.23 


151.22 


AMI-LA 


15.4 


29.31 ±0.73 


152.27 


AMI-LA 


15.4 


28. 10 ±1.30 



TABLE 1 — Continued 



St" 


Facility 


Frequency 


nux Density 






(GHz) 




153.46 


AMI-LA 


15.4 


26.71 ±0.38 


155.45 


AMI-LA 


15.4 


27.80 ± 0.45 


156.35 


AMI-LA 


15.4 


29.92 ±0.89 


157.45 


AMI-LA 


15.4 


27.77 ±0.31 


158.29 


AMI-LA 


15.4 


27.72 ±0.83 


159.28 


AMI-LA 


15.4 


27.98 ±0.59 


160.27 


AMI-LA 


15.4 


28.93 ±0.44 


161.20 


AMI-LA 


15.4 


26.94 ±0.84 


162.21 


AMI-LA 


15.4 


26.81 ±1.21 


163.28 


AMI-LA 


15.4 


28.82 ±0.68 


166.00 


AMI-LA 


15.4 


26.92 ±0.58 


168.13 


AMI-LA 


15.4 


25.90 ±0.64 


169.13 


AMI-LA 


15.4 


25.47 ±0.83 


172.25 


AMI-LA 


15.4 


26. 12 ±0.48 


174.33 


AMI-LA 


15.4 


26.12±0.15 


177.32 


AMI-LA 


15.4 


26.31 ±0.37 


181.30 


AMI-LA 


15.4 


24.37 ±0.85 


184.03 


AMI-LA 


15.4 


23.06 ±0.77 


187.19 


AMI-LA 


15.4 


24.24 ±0.43 


187.91 


AMI-LA 


15.4 


24.28 ± 0.45 


189.37 


AMI-LA 


15.4 


24.60 ± 0.59 


193.19 


AMI-LA 


15.4 


26.15 ±0.75 


194.01 


AMI-LA 


15.4 


24.53 ±0.34 


197.02 


AMI-LA 


15.4 


26.72 ±0.52 


200.01 


AMI-LA 


15.4 


25.90 ±0.15 


203.01 


AMI-LA 


15.4 


24. 12 ±0.41 


207.08 


AMI-LA 


15.4 


25.06 ±0.45 


210.19 


AMI-LA 


15.4 


24.24 ±0.31 


214.17 


AMI-LA 


15.4 


22.64 ±0.38 


217.16 


AMI-LA 


15.4 


24. 16 ±0.25 


220.20 


AMI-LA 


15.4 


23.74 ±0.05 


221.22 


AMI-LA 


15.4 


24.33 ± 0.49 


225.91 


AMI-LA 


15.4 


23.80 ±0.14 


232.87 


AMI-LA 


15.4 


23.72 ±0.16 


235.01 


AMI-LA 


15.4 


21.64±0.16 


237.85 


AMI-LA 


15.4 


20.86 ±0.75 


4.79 


EVLA 


19.1 


2.12±0.02 


6.75 


EVLA 


19.1 


4.36 ±0.05 


7.77 


EVLA 


19.1 


5.25 ±0.03 


8.87 


EVLA 


19.1 


6.38 ±0.06 


9.78 


EVLA 


19.1 


5.42 ±0.03 


21.89 


EVLA 


19.1 


12.01 ±0.03 


31.74 


EVLA 


19.1 


13.50 ±0.05 


35.86 


EVLA 


19.1 


13.97 ±0.05 


50.65 


EVLA 


19.1 


17.11 ±0.06 


67.61 


EVLA 


19.1 


23.03 ±0.06 


94.64 


EVLA 


19.1 


31.36 ±0.07 


111.62 


EVLA 


19.1 


30.21 ±0.10 


127.83 


EVLA 


19.1 


29.75 ±0.22 


142.62 


EVLA 


19.1 


29.57 ±0.13 


159.77 


EVLA 


19.1 


26. 10 ±0.26 


177.50 


EVLA 


19.1 


24.24 ±0.16 


198.22 


EVLA 


19.1 


23.02 ±0.12 


219.22 


EVLA 


19.1 


23.15 ±0.07 


4.79 


EVLA 


24.4 


3.01 ±0.03 


6.75 


EVLA 


24.4 


5.58 ±0.06 


7.77 


EVLA 


24.4 


6.70 ±0.03 


8.87 


EVLA 


24.4 


7.88±0.12 


9.78 


EVLA 


24.4 


6.84 ±0.03 


21.89 


EVLA 


24.4 


12.69 ±0.02 


31.74 


EVLA 


24.4 


13.80±0.06 


35.86 


EVLA 


24.4 


14.95 ±0.05 


50.65 


EVLA 


24.4 


18.30±0.06 


67.61 


EVLA 


24.4 


25.62 ±0.06 


94.64 


EVLA 


24.4 


30.67 ± 0.07 


111.62 


EVLA 


24.4 


28.20 ±0.13 


127.83 


EVLA 


24.4 


28.29 ±0.33 


142.62 


EVLA 


24.4 


24.73 ±0.16 


159.77 


EVLA 


24.4 


23.83 ±0.40 


177.50 


EVLA 


24.4 


20.40 ±0.19 


198.22 


EVLA 


24.4 


19.88±0.15 


219.22 


EVLA 


24.4 


21.40 ±0.08 


5.75 


EVLA 


43.6 


8.11 ±0.16 


6.75 


EVLA 


43.6 


7.70 ±0.14 


8.87 


EVLA 


43.6 


9.62 ±0.14 



TABLE 1 — Continued 



St"" 


Facility 


Frequency 


Flux Density 


(d) 




(GHz) 


(mJy) 


111.62 


EVLA 


43.6 


22.09 ± 0.65 


127.83 


EVLA 


43.6 


20.84 ± 0.48 


159.77 


EVLA 


43.6 


15.43 ±0.63 


177.50 


EVLA 


43.6 


LS.73±0.38 


198.22 


EVLA 


43.6 


1L70±0.26 


219.22 


EVLA 


43.6 


15.32±0.15 


4.90 


CARMA 


87 


15.66±0.51 


8.19 


CARMA 


87 


20. 16 ±0.96 


9.14 


CARMA 


87 


21.67 ±0.33 


10.23 


CARMA 


87 


14.70 ±1.28 


12.14 


CARMA 


87 


17.67 ±0.94 


1 7 Aft 




o / 


IZ. 1 J It 1 . JZ 


90 1 1 


PARMA 


o / 


11 70-1-1 74 


ZJ. IZ 


PARMA 


o / 


1D.7J HI J.^o 


QQ 7^ 


PARMA 


R7 


1 Q 70 -L n rt*? 


1 J L.DZ 


P AT?MA 


S7 


1 1 S -t- ft '^'^ 
1 J. 1 o III yf.JD 






o / 


i 1 .oo m U.J J 


1 / J. jD 


PARMA 


R7 


y.jy m u. J 1 


1 Q 9^ 




zuu 


1 A 1 n 4- 1 

l^.lUlt I.JU 


24.32 


SMA 


200 


10.70 ±1.00 


10.30 


SMA 


230 


14.90 ±1.50 


11.13 


SMA 


230 


11.70±1.40 


17.23 


SMA 


230 


13.30±1.50 


18.25 


SMA 


230 


9.90 ±1.40 


20.24 


SMA 


230 


8.20 ±1.40 


21.25 


SMA 


230 


8.30 ±2.20 


125.05 


SMA 


230 


6. 10 ±0.65 


5.13 


SMA 


345 


35. 10 ±0.80 



Note. — 



" AH values of St are relative to the initial 7-ray detection: 2011 March 25.5 UT. 



TABLE 2 

Results of Broad-band Spectral Energy Distribution Fits 



St 


log(l^a) 


log(f„) 




log(rig) 


log(rsh) 


log(r;) 


log(Lj iso,48) 


log(ni8) 


log(/7cNM) 


log(«j) 


(d) 


(Hz) 


(Hz) 


(mJv) 


(cm) 






Ter!? s ' ^ 


(cm^^) 


(cm"') 


(cm"') 


5 


11.01 


11.74 


1.47 


-0.79 


0.65 


0.78 


-1.33 


0.28 


1.82 


3.53 


10 


10.11 


11.41 


0.96 


-0.26 


0.61 


0.74 


-1.39 


-0.28 


0.25 


2.39 


15 


10.02 


11.21 


0.98 


-0.13 


0.55 


0.65 


-1.33 


-0.22 


0.05 


2.19 


22 


9.96 


10.99 


0.97 


-0.04 


0.51 


0.60 


-1.25 


-0.15 


-0.07 


2.08 


36 


9.95 


10.78 


1.11 


0.08 


0.47 


0.54 


-1.05 


0.01 


-0.17 


2.03 


51 


9.96 


10.62 


1.21 


0.16 


0.43 


0.50 


-0.92 


0.13 


-0.18 


2.01 


68 


9.99 


10.53 


1.39 


0.24 


0.41 


0.47 


-0.72 


0.28 


-0.20 


2.04 


97 


9.95 


10.39 


1.52 


0.36 


0.40 


0.45 


-0.53 


0.39 


-0.34 


1.99 


126 


9.97 


10.26 


1.58 


0.41 


0.36 


0.41 


-0.46 


0.49 


-0.33 


1.97 


161 


9.82 


10.13 


1.51 


0.52 


0.36 


0.40 


-0.39 


0.44 


-0.60 


1.82 


197 


9.83 


10.04 


1.56 


0.56 


0.34 


0.38 


-0.32 


0.52 


-0.60 


1.81 


216 


9.90 


9.99 


1.63 


0.55 


0.32 


0.35 


-0.29 


0.60 


-0.50 


1.85 



Note. — Infen'ed parameters of the relativistic outflow and eiiviroimient of Sw 1644+57 from model fits of individual multi-frequency epochs. The model 
is described in 
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10 10 10 10 10 10 

Time (d) Time (d) 



Fig. 1 . — Radio light curves of Sw 1644+57 extendingto (5f Ri 216 d. The data at (5f Rd 5 -22 d were previously presented in ZBSl 1. The solid lines are models 
based on independent fits of broad-band SEDs (Figure [2) using the model described in 0(see also MGMll). The dashed lines are the predicted light curves 
from MGMll, which assumed a constant energy and a steady density profile of p oc f^. The secondary maximum in the millimeter band and the continued 
increase in brightness to a peak time of St > 100-200 d in the centimeter bands require an increase in the energy and density relative to the initial evolution. 
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Fig. 2. — Radio spectral energy distributions of Sw 1644+57 at 5? Ri 5 — 216 d. The data at 5; 5 -22 d were previously presented in ZBSl 1. The solid lines 
are fits based on the model described in ^(see also MGMll). In each epoch we fit forLj iso and nig with fixed values of = eg = 0.1 and p = 2.5. The dashed 
gray fines mark the peak flux density and peak frequency at <5f = 10 d to help track the evolution of the spectrum as a function of time. The dashed red lines mark 
the expected SEDs based on the evolution at St ^ 5 — 22 d. 
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Fig. 3. — Model and extracted parameters for each broad-band SED shown in Figure[2] Shown are the time evolution of the synchrotron parameters (Va, Um, 
and /Vo), Lj isoAS' f^is. ^sh, and Tj. The substantial increase in energy and density is clearly seen. In addition, we find r oc r"^, a steeper increase than expected 
in a simple Wind model with constant energy 
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Time (d) 

Fig. 4. — Temporal evolution of the synchrotron frequencies v,„ (top) and Va (bottom) relative to the expected evolution in a simple model with a constant 
energy and a Wind profile. The shallower dechne of oc ii^'^, with a particularly shallow evolution at 6t ~ 30—100 d, is indicative of a continuous increase in 
energy. Similarly, the shallower decline of Va, followed by a rapid increase, is indicative of a density profile of p oc r~'^ and a flattening at i5f « 30- 100 d. 
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Time (d) 

Fig. 5. — Temporal evolution of the integrated luminosity (or alternatively £j,iso = ij.isoO; black circles) in comparison to the X-ray luminosity (gray dots) 
as parametrized with a simple luminosity evolution (gray line). The red curve is the integrated luminosity derived from the simple model. The observed X-ray 
luminosity indicates that the fraction of total energy emitted in X-rays is comparable to the energy in the relativistic outflow (i.e., ex ~ 0.5). The large increase 
in energy infen'ed from the radio observations cannot be explained by injection from a L oc r^/^ tail. 
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Fig. 6. — Radial density profile in the inner parsec around Sw 1644+57 as inferred from the radio observations (black circles). The overall profile follows 
p oc r"^!"^, with a significant flattening at r 0.4-0.6 pc. Following the flattening, the profile appears to recover to r"^!- by about 1 pc. Also shown is the density 
infen'ed from X-ray observations of the Galactic center (gray squares; Baganoff et al. 2003), which is about a factor of 30 times larger at si 0.5 pc. The dashed 
line is a scaled-down model of the Galactic center density profile assuming gas feeding from massive stars in which the bulk of the gas is thermally expelled in a 
wind JOuataert.i2004) . In this model the inner profile (< 0.2 pc) is oc r"'^!'^, while the outer profile (> 0.4 pc) has a Wind (oc r"^) profile. 
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Time (d) 



Fig. 7. — Predicted optical (r-band; blue) and near-infrared (A'-band; red) light curves using the results of the radio modeling. The upper limits in r-band and 
detections in S'-band are from Levan et al. (2011). Since the A'-band fluxes are the total for Sw 1644+57 and its host galaxy we have subtracted an estimated 
host contribution of about 20 fi]y {K 20.6 AB mag). The solid lines are models without a cooling break between the radio and optical/near-IR, which clearly 
over-estimate the AT-band flux. The dashed lines include a cooling break at Uc ~ lO" Hz, and the dotted lines add host galaxy extinction of Av ~ 3.5 mag to 
account for the optical non-detections. The combination of a coohng break and extinction provides an excellent fit to the near-IR evolution. 
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Fig. 8. — Predicted evolution of the radio liglit curves at 6 GHz (blue) and 22 GHz (red) assuming a radial density profile of p oc r at r > 1 pc and three 
values for the maximum integrated beaming-corrected energy (solid: Ej = 10^- erg; dashed: Ej = 3 X 10^^ erg; dotted: Ej = 10^^ erg). The thin horizontal 
hnes mark the 5a sensitivity of the EVLA, and indicate that the radio emission from Sw 1644+57 should be detectable for decades (and perhaps centuries) at 
centimeter wavelengths. 
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Fig. 9. — Predicted evolution of tlie jet radius assuming a radial density profile of p oc r '-^ at r > 1 pc and three values for the maximum beaming energy 
(solid: Ej = 10^- erg; dashed: Ej = 3 X 10^~ erg; dotted: Ej = lO'^ erg). The thin horizontal line marks the resolution of VLBI for ajet opening angle of dj = 0.1. 
The source should become resolvable at (5r ~ 6 yr, with an expected 22 GHz flux density of about 2 mjy (Figure[8). If the jet instead begins to undergo significant 
spreading it may become resolvable at ~ 1 yr when the 22 GHz flux density is still ~ 10 mJy. 



